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physical disks. Although the hydro dynamic stability of disk-like flows has been tested in exper- 
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Abstract 

^^. Turbulent transport of angular momentum is a necessary process to explain accretion in astro- 

i 

C \ iments, results are contradictory and suggest either laminar or turbulent flow. Direct numerical 

O ■ simulations reported here show that currently investigated laboratory flows are hydrodynamically 

^. unstable and become turbulent at low Reynolds numbers. The underlying instabilities stem from 

Q^\ the axial boundary conditions, affect the flow globally and enhance angular momentum transport. 
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Accretion in astrophysical disks requires the flow of mass towards a central gravitating 
body. The ensuing loss of momentum must be balanced by outward angular momentum 
transfer among gas particles [1]. If the motion of orbiting gas was laminar molecular trans- 
port would be orders of magnitude too slow for accretion to take place, and so considering 
a turbulent viscosity becomes necessary [2], However, in Keplerian disks the gas rotates as 
Q oc r~^/^ and laminar motion is linearly stable according to the Rayleigh criterion. Al- 
though axial magnetic flelds can drive turbulence via the magnetorotational instability [3|, 
it is not clear whether this operates in weakly ionized disks. On the other hand, it is well 
known that linearly stable shear flows (such as pipe flow) can become turbulent due to flnite 
amplitude disturbances. Whether Keplerian flows are susceptible to such a transition sce- 
nario or remain stable despite the large Reynolds numbers (i?e), is a topic of great interest 
and the source of much controversy ^] . 

The stability of disk-like flows is typically probed in laboratory experiments of fluid 
between two concentric and independently rotating cylinders, Taylor-Couette flow (TCF). 
In the inflnite-cylinder idealization, the Navier-Stokes equations admit a pure rotary solution 

commonly referred to as Couette flow. Here ri and r2 are the radii of the inner and outer 
cylinders and ^i and 0.2 their angular velocities. When {rxjr^)^ < ^^2/^1 < 1 the angular 
velocity decreases outward but the angular momentum increases (quasi-Keplerian flows). 
Accretion disks are stratifled in the axial direction and are thus best modelled considerine; 
an abounded d„.a.„. thereby av„.d.g .t.fic.al bou.da,, c„„d,t.„„s g Most expen- 
ments and simulations focus, however, on the physics of the disk's midplane and neglect 
stratiflcation. Under this assumption simulations typically employ axially periodic bound- 
ary conditions, whereas in experiments cylinders have a flnite-length h. Hence the degree 
to which ([T]) may be approximated is compromised by the axial boundary conditions and 
length-to-gap aspect-ratio F = h/{r2 — ri). In particular, solid axial boundaries result in a 
basic state with nonzero radial and axial velocity components (Ekman flow). Hence, pro- 
ducing Couette-like proflles in experiments poses an extraordinary challenge, which may be 
addressed by considering very tall cylinders or splitting the endwalls in several rings that 
rotate at independent angular speeds M. The latter strategy has been implemented in the 
Princeton Taylor-Couette experiment [Z], which has a short aspect-ratio F = 2.104 but 



whose endwalls are split into two independently rotating rings. From simultaneous Laser 
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Doppler Velocity measurements of azimuthal and radial velocity components, Ji et al. J8(] 
suggest that quasi-Keplerian flows at Reynolds numbers up to millions are essentially lam- 
inar. From this, they conclude that purely hydrodynamic mechanisms cannot transport 
angular momentum at the rates required for accretion to occur in disks. 

This conclusion has been recently challenged in a new experimental study by Paoletti and 
Lathrop [9] , who report from direct torque measurements at the inner cylinder that Keplerian 
flows at i?e > 10^ are fully turbulent. When extrapolated to astrophysical disks, their 
results indicate that transport occurs at accretion relevant rates, in agreement with previous 
estimations J10l|. Despite having tall cylinders F = 11.47, the Maryland experiment has solid 
endwalls that are attached to the outer cylinder and hence cannot be rotated independently. 
Although these boundary conditions are known to generate vigorous Ekman vortices and 
greatly increase the exerted torque, their contribution is discarded by dividing the inner 
cylinder into three sections and measuring torque only in the central one. Despite efforts in 
the Princeton and Maryland experimental setups to mitigate endwall effects it is, however, 
unclear whether their results can be used to infer the stability of flows in astrophysical disks 
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J]. In this Letter it is shown that current laboratory experiments of quasi-Keplerian flows 
become turbulent already at Re = O(IO^) due to hydrodynamic instabilities stemming from 
the axial boundary conditions. Moreover, it is found that turbulence flUs the entire flow 
domain and as a result the momentum transfer is globally enhanced. 

Here direct numerical simulations of flows with the precise geometry and boundary con- 
ditions of the Princeton and Maryland experiments were performed. The Navier-Stokes 
equations for an incompressible Newtonian fluid of velocity v 

dtv + (v • V)v = -Vp + Av, V • V = , (2) 

were rendered dimensionless by scaling lengths and time with the gap-width d = r2 — ri 
and viscous time cP/u^ where u is the kinematic viscosity of the fluid. The solution of ([2j) 
was formulated in primitive variables in cylindrical coordinates (r, ^, z) and a second-order 



time-splitting method with consistent boundary conditions for the pressure was used |ll| . 
The spatial discretization consists of Chebyshev collocation in (r, z) and a Galerkin- Fourier 
expansion in 9. The code converges spectrally in the three directions [12| and was validated 
against a Legendre-Fourier-Galerkin code p^]. Here the resolution was chosen to ensure 



that computed torque values were accurate to at least 1%. 

The geometry of the Taylor-Couette system is specified by the radii-ratio r] = ri/r2 and 
the length-to-gap aspect-ratio F. The dimensionless boundary conditions at the cylinders 
read {vr^vo^Vz)[ri^2^0^ z] = (0, i?ei,2, 0), where Rei = driVti/v {Re2 = dr2^2/i^) is the inner 
(outer) cylinder Reynolds number. Because of differential rotation the angular velocity 
changes abruptly at adjacent rotating boundaries. In the Princeton experiment, the endwall 
is split at mid-radius r^ = (ri+r2)/2 into two independently rotating rings. Hence, there are 
four independent angular speeds; Qi and f^2, for inner and outer cylinder, and Qs and 1^4, for 
inner and outer rings. To preserve spectral convergence discontinuities in angular velocity 
were regularized, yielding the following boundary condition at the endwalls z = ±r/2 

n{r)={n^-ns)exp[-{r-n)/e] + 
(f^2-^4)exp[-(r2-r)/e] + 

^ + ^ tanh[(r - r^)/e], 

with 6 G [5 X 10~^, 10~^] (see Ref. 14). The boundary condition modeling the Princeton 
experiment [8] is shown as circles in Fig. [It. Due to the sharp gradient dQ/dr\r^ and 
the clustering of Chebyshev points close to the boundaries, a large number of radial points 
{rir = 351) was required to accurately simulate the split endwall. In the axial and azimuthal 
directions up to n^ = 281 Chebyshev points and uq = 256 Fourier modes were used. The 
Maryland experiment has a single solid ring attached to the outer cylinder (Q^ = Q^ = Q2) 
and there is only a strong gradient at ri (see crosses in Fig. [lb). Here up to n^ = 601, 
no = 384 and n^ = 61 were used. 

To put TCF in the wider context of rotating shear fiows it is useful to define a shear 
Reynolds number i?e = 2/(1 + r])\Re2r] — Rei\ and a rotation number i?^ = (1 — r]){Rei + 
Re2)/{Re2Ti — i?ei), which measure the ratio of shear to viscous forces and the ratio of 
mean rotation to shear [151] , respectively. Here the sign of R^ distinguishes between cyclonic 
(i?^ > 0) and anticyclonic fiows (i?^ < 0), with —2 < i?^ < — 1 corresponding to quasi- 
Keplerian rotation. In the experiments of Ji et al. [8] i?^ = —1.038 and tj = 0.3478, and the 
same values were used here, whereas Lathrop and Paoletti [9] have systematically studied 
both cyclonic and anticyclonic regimes at ry = 0.7245. Here, r] = 0.7245 and the value Rq = 
— 1.047 was chosen (corresponding to their Rossby number Ro = Rei/{r]Re2) — 1 = 0.85). 

The endwall infiuence in the Maryland experiment is illustrated in Fig. [1^, showing the 
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FIG. 1. (color online) (a) Angular velocity at the endwalls for the Princeton (circles) and Mary- 
land (crosses) experiments, and Couette flow (dashed lines), (b-c) Steady basic states at Re = 320 
(Maryland) and Re = 772 (Princeton). White (black) corresponds to maximum (minimum) veloc- 
ity and the radial direction is horizontal with left (right) corresponding to inner (outer) cylinder. 

velocity field of the steady basic state at Re = 320. At the endwalls there is a strong 
negative radial velocity infiow, which generates axial velocities pointing towards mid-height 
along the inner cylinder and result in an axially dependent azimuthal velocity. Figure [lb 
shows the basic state of the Princeton experiment at Re = 772. Because of the small 
aspect-ratio the meridional circulation generates a strong radial outward fiow at mid-height 
that increases the outward transport of azimuthal velocity. Were the anerular speeds of the 

nn 

endwall rings selected according to the ideal Couette profile [y, 116|] instead of the values 
used in experiments [8] and reproduced here, profiles with weaker meridional circulation 
and hence closer to Couette fiow could be obtained. This approach was used in previous 
numerical simulations of the Princeton configuration \X^ . 

The visualizations of Fig. [T^-c hint at the difficulty of realizing quasi-Keplerian profiles 
in a laboratory experiment even at very low Reynolds numbers. The endwall boundary 
conditions change the velocity field not just locally but globally across the domain. In fact. 




FIG. 2. (color online) Three-dimensional view of isosurfaces of negative (red) and positive (yellow) 
radial velocity, (a) tti = 2 rotating wave at Re = 1332 and (b) turbulent flow at Re = 5328 from 
simulations of the Maryland experiment, (c) Modulated rotating wave with m — 1 and m — 2 si 
Re — 1545 and (d) turbulent flow at Re — 6437 from simulations of the Princeton experiment. 



at slightly higher Reynolds numbers the flow becomes three-dimensional and time-dependent 
via supercritical Hopf bifurcations. In the Maryland configuration instability occurs first at 
Rcc = 352 to prograde rotating waves with azimuthal wavenumber m = 5 and localized at 
the endwalls. Beyond Rcc multiplicity of states, with different symmetries and wavenumber 
m G [2, 5], was found, whereas for Re > 1330 only a global m = 2 mode remained stable and 
was obtained regardless of initial conditions (see Fig. [Sb). Further increasing the Reynolds 
number led to modulated waves and a quick transition to temporal chaos at about Re c^ 1600. 
Subsequently, spatial periodicity was lost and the spectra broadens as the fiow became 
gradually turbulent; Fig. [2t shows a fiow snapshot at Re = 5328. This transition picture 
is also representative of the simulations of the Princeton experiment. Here the basic steady 
state becomes unstable at Re = 1448 almost simultaneously to m = 1 and m = 2 rotating 
waves, which were found to coexist in space and time (see Fig. [2b). By no means is this 
situation generic: changing the relative rotation of the cylinders one of m = 1, 2 was found 
to bifurcate first. Further increasing the Reynolds number led to very complex and strongly 
three-dimensional fiow as shown in Fig. [2]i. 

The stability of quasi-Keplerian TCF with endwalls attached to the outer cylinder and 
T] = 0.7245 is shown in Fig. [3h. The minimum critical Reynolds number is attained at the 
Rayleigh line {Rq = — 1) and increases as diflFerential rotation decreases, but with Rcc < 10^ 
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FIG. 3. (color online) Stability curves of quasi-Keplerian Taylor-Couette flows with endwalls 
attached to the outer cylinder: (a) 77 = 0.7245 and F as in the legend, (6) R^ — —1.1 and 77 as in 
the legend. 



across the whole quasi-Keplerian regime. Close to the Rayleigh line modes localized to the 
endwalls bifurcate first, whereas for i?^ < —1.1 global instability modes (as in Fig. [2h) 



18 and the latter 



dominate. The former are similar to those observed experimentally in Ref. 

are similar to those reported by Avila et al. il2|, who studied global boundary layer effects 

on fiows between exactly co-rotating cylinders and stationary endwalls. Figure [St further 

shows that endwall instabilities depend weakly on geometry and hence generically govern 

the dynamics of quasi-Keplerian TCF. It is worth noting that endwall instabilities persist 

beyond the Rayleigh line and coexist with Taylor vortices close to the onset of centrifugal 

instability. 

The onset of hydrodynamic instability and transition to turbulence are expected to rad- 
ically change the radial transport of azimuthal momentum. The solid lines in Fig. Hb show 
normalized average azimuthal velocity profiles {ve)e,t/{'^i^i) at mid-height for simulations 
of the Maryland experiment at Re = 5328 (black solid line) and Princeton experiment at 
Re = 6437 (gray solid line, orange online). At the inner cylinder the profiles are steeper 
than laminar Couette fiow (dashed lines), implying larger torques on the cylinder surface. It 
is worth noting that in TCF between infinite cylinders the transverse current of azimuthal 
motion J^ = r^[{vr^)o^z,t — ^dr{^)o^z,t\ is a conserved quantity [19], and as a consequence 
the dimensionless torque G = iy~'^J^ is the same at inner and outer cylinder. This does not 
hold, however, for fiows confined by no-slip axial boundaries. Torque profiles along the inner 
cylinder, normalized with the laminar Couette torque, are shown in Figure Hf). Because of 
the sharp change in Q occurring across a small gap between inner cylinder and endwalls 
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FIG. 4. (color online) (a) Average azimuthal velocity at mid-height for simulations of the Maryland 
{Re = 5328, black solid line) and Princeton (Re = 6437, orange solid line) experiments. The dashed 
lines are laminar Couette flow, (b) Torque Nu = G/Giam along the inner cylinder, solid lines as 
in (a), (c) Reynolds number dependence of the torque on the central section of the inner cylinder 
(0.6 > |2z/r| and 0.4 > |22:/r| for simulations of the Maryland and Princeton experiments). 



(see Fig. [Tb), the torque required to rotate the inner cylinder faster than the endwall is 
very large. This increase in local torque as the endwalls 2z/F = ±1 are approached can be 
seen in Figure H^. Although the direct contribution of the endwall is largely avoided by the 
measurement technique in the experiments, the torque in the central section remains well 
above laminar because of turbulent fluctuations. 

The simulations of the Maryland experiment show a clear change in the torque behavior 
at about Re = 3000 (see black curve in Fig. \^). This is related to the appearance of the two 
torque peaks at 2z/r ^ ±0.5 in Fig. Hf) and is caused by the onset of small-scale vortices 
which have opposite spiral orientation to the structure of the primary rotating wave (see 
FigEJ)). In both experiments the torque has already doubled the laminar value at Re ^ 6000 
due to the endwall driven instabilities. 

In spite of the disparity in Reynolds numbers, it is tempting to compare these numerical 



results to experimental observations. Van Gils et al. J20| have measured torque at Re > 10^ 
and have observed an effective universal scaling law Nu = G/Giam = a{R^)Re^'^^ that 
'lolds throughout the linearly unstable regime of TCF. Surprisingly, Paoletti and Lathrop 
9| (see also Ref. |2l|) have demonstrated that this law applies to linearly stable cyclonic and 
anti-cyclonic Rayleigh-stable regimes as well, that is including quasi-Keplerian rotation. It 
is then natural to ask how this universal behavior connects to the complex flows uncovered 
in this work. It is speculated here that a transition between end wall- driven turbulence to 
the universal Nu oc Re^'^^ scaling reported in experiments may take place at intermediate 
Re. Interestingly, in the case of stationary outer cylinder a cross-over marking the transi- 
tion from centrifugally to shear-driven turbulence at Re c^ 13000 was reported J22l]. If an 
analogous cross-over was found in quasi-Keplerian flows and shown to be independent of 
aspect-ratio and endwall boundary condition, a strong case for the existence of hydrody- 
namic turbulence in astrophysical disks would be made. In fact, ingredients of shear-driven 
turbulence such as transient growth of disturbances are found also in quasi-Keplerian flows, 
although signiflcantly only at Re = O(IO^) [23]. On the other hand, it would be interesting 
to investigate the connection between the complex flows found here and the quiescent flows 
reported by Ji et alfs] at large Re. 

In conclusion, current laboratory experiments designed to approximate flow proflles ex- 
pected from accretion disks become turbulent at moderate Reynolds number due to imposed 
boundary conditions. Although these instabilities are generic and hence cannot possibly be 
avoided, universal scaling suggests that shear mechanisms might overwhelm endwall eflFects 
at large Reynolds number. In order to probe this hypothesis new experiments with variable 
aspect ratio and diflFerent axial boundary conditions should be conducted. These would 
provide great indsight on the physical mechanisms of rotating shear flows and might shed 
light on the origin of turbulence in astrophysical disks. 
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